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ABSTRACT 

Recent inspections of local available data suggest that the almost linear relation between the 
stellar mass of spheroids (Msph) and the mass of the super massive Black Holes (BHs) resid¬ 
ing at their centres, shows a break below Mgph ^ 10^° M©, with a steeper, about quadratic re¬ 
lation at smaller masses. We investigate the physical mechanisms responsible for the change in 
slope of this relation, by comparing data with the results of the semi-analytic model of galaxy 
formation MORGANA, which already predicted such a break in its original formulation. We 
find that the change of slope is mostly induced by effective stellar feedback in star-forming 
bulges. The shape of the relation is instead quite insensitive to other physical mechanisms 
connected to BH accretion such as disc instabilities, galaxy mergers. Active Galactic Nucleus 
(AGN) feedback, or even the exact modelling of accretion onto the BH, direct or through a 
reservoir of low angular momentum gas. Our results support a scenario where most stars form 
in the disc component of galaxies and are carried to bulges through mergers and disc insta¬ 
bilities, while accretion onto BHs is connected to star formation in the spheroidal component. 
Therefore, a model of stellar feedback that produces stronger outflows in star-forming bulges 
than in discs will naturally produce a break in the scaling relation. Our results point to a form 
of co-evolution especially at lower masses, below the putative break, mainly driven by stellar 
feedback rather than AGN feedback. 
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1 INTRODUCTION 


Since 1998, observations have shown the existence of more or less 
tight correlations between the mass of the super-massive Black 
Holes (BHs), hosted at the centre of local galaxies, and sev- 
eral properties of their spheroidal component like luminosity (e.g . 
iMagorrian et alJll 998l:lMarconi & HuiHl2003l: lLasker et'^l2014 ). 


stellar mass (e.g. Haring & Rixll2004 ; SanLgt_alJ|20ilJ), i^tral ve 

locity dispersion (e.g. llVterritt & Ferraresel 1200 ll: iTremaine et al.l 


locity dispersion (e.g. [ 
2003: IShankail 120091) 


Graham et al.ll20oi : Savorgnan et alj2013h . 


or radial concentrate of stars (e.g. 


The discovery of these scaling relations has been the start¬ 
ing point for a re-evaluation of the role of accretion of gas into 
BHs and their energetic feedback due to Active Galactic Nu¬ 
cleus (AGN) activity, in the context of galaxy evolution. Start- 
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l2012l : lHirschmann et alj2()T^ : lMenci et al.l2014h have investigated 
the physical mechanisms responsible for these tight relations. Most 
of these studies have been developed in the framework of “co¬ 
evolution” of BH growth and host galaxies, where either the modu¬ 
lation of accretion or AG N feedback are responsible for the co rrela- 
tions. Other authors (e.g. lPen 32007l : ljahnke & Maccicll201 ll) have 
instead pointed out that a linear relation between bulge and BH 
masses is not necessarily a sign of co-evolution, but is also com¬ 
parable with a simple “co-habitation” model where correlations are 
mainly driven by galaxy mergers, via the central limit theorem. 

In the last years, new observations have revisited and some¬ 
times modified these scaling relations. In particular, recent up¬ 
dates have been published pointin g to higher normalizations in the 
BH m ass-bulge stellar mass (e.g.. lGrahamll2013 : lKormendv & Hoi 
2013h and BH mass-central velocity dispersion relations (e.g., 
GrahamiboiSh . This in turn could imply significantly larger to¬ 
tal m ass densities for the super-massive BHs in the loc al Universe 
(e.g.. lNo\^l2013l : IShankaill2013l : IComastri et alj|201^ . These re¬ 
visions are due to a number of improvements both in the obser¬ 
vational data, thanks to adaptive optics and integral-field spec- 
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troscopy, and in the modeling of star kinematics. Bulge masses are 
intrinsically more difficult to determine, with respect to other bulge 
properties like luminosity or velocity dispersion, due to the uncer¬ 
tainties in mass-to-light ratio (which stem from the uncertainties 
in the initial stellar mass function, the modelling of stellar evolu¬ 
tion and the age-metallicity degeneracy). A breakthrough in this 
respect came from the s tudy of galaxy kinematics in SAUROIsQ 
and ATLAS-3D surveys jCaDpellari et alJl2006Ll2013h . which pro¬ 
vided calibrated relations between mass-to-light ratios and, e.g., ve¬ 
locity dispersion for a representative sample of local galaxies. All 
these developments led to a significant recali bration of the ov er- 
all MBn/Afsph ratio, from ~ 0.1 perce nt (e.g. lSani et alJl20TIh to 
0.49 percent (e.g iGraham & Sc^l2013h . 

Moreover, some evidence has also been accumulating sug¬ 
gesting a scenario where different galaxy populations may follow 
individual (and different) relations. For example, pseudo-bulges, 
usually defined as spheroids characterized by disk-like exponential 
profiles and/or rotational kinematics, were reported to systemati¬ 
cally lie below the main relation defin ed by classical bulges and 
ellipticals (see, e.g ., the discussion in Gr^aml boOSl : IHii l2008l : 
IShankar et al.ll20la : iKormendv & Hdl2013l) . but recent studies do 
not confirm these findings (Savorgnan et ah, in preparation). The 
power-law, almost linear relation (in logarithmic space) between 
BIT mass and either bulge stellar mass or central velocity dis¬ 
persion has been often used as a constraint for models of the 
joint evolution of galaxies and AGNs. Indeed, although theoret¬ 
ical models may trace diffe rent evolutionary path s for different 
galaxy populations (see, e.g.. lLamastra et ahllTOlOll . they typically 
predict a local Mhh-Msph relation compatible with a single power- 
law down to small mass scales that are barely probed by observa¬ 


tions (e.g.. iMa rulli et al 


I 2 OI 2 I : iHirschmann et al 


200 s; iMenci et alj|2008| : iFanidakis et al.l 


2012n . A remarkable exception is pro¬ 


vided by the MOde l for Rise of GAl axies aNd AGNs (MO RGANA, 
[Monaco et al.l200^ . and presented in lFontan^ et al.l ^20061) . In that 
paper, the predicted Mhh-Msph relation had a change of slope at 
Mbh ~ 10® Mq, with a steeper relation at smaller masses. This 
break had no observational support, but was not incompatible with 
the sparse data available at that time. 

In what follows, whenever we refer to the Black Hole - Bulge 
relation (Mhh-Msph relation), we will specifically focus on the 
relation between BH mass (Mbh) and the stellar mass of the 
spheroidal/bulge component of the host galaxy (Msph). 

The overall shape of the Mhh- Msph relation over a wide range 
in Msph has been recently revised bv iGraha m & Scot^(|2015l), wh o 
studied the Mhh-Msph relation using data from Scott et alj( l2013b . 
and included AGN data for which Mgph had been derived for the 
first time. These authors found that the relation is linear for BH 
masses above Mbh ^ 2 x 10® Mq, while it significantly steep¬ 
ens towards a quadratic relation below this threshold. They argued 
that previous studies missed the “break” in the Mhh-Msph relation 
due to an insufficient sampling at the low-mass end, mainly below 
Mbh lO^M©. It must be however stressed that consensus on the 
statistical r elevance of the break in Mhh-Msph relation has still to 
be reached. iGraham & Scotll ( 1201 Sh further revealed that the most 
massive galaxies defining the linear part of the Mhh-Msph relation 
are mostly described by a core-Sersic stellar profile, i.e a profile 
with a deficit of light in central regions with respect to the extrap¬ 
olation of their outer profile, while the less massive ones are better 
defined by a single Sersic profile. 


^ Spectrographic Area Unit for Research on Optical Nebulae 


The existence of a quadratic, if not steeper, Mhh-Msph rela¬ 
tion has key implications for the modeling of the joint evolution 
of BHs and their host galaxies. In particular, it implies that the 
Mbh /Msph ratio is not constant, but it rather depends on the final 
amount o f stellar mass in t he ho st galaxy, and/or its morpholog¬ 
ical type. iGraham & ScottI ( 1201 5h proposed a physical interpreta¬ 
tion of their results suggesting that the initial dissipative processes 
controlled by the gas-rich, initial phases of galaxy-BH formation, 
establish a quadratic relation. Later, mostly dry (gas-poor) mergers 
are responsible for the gradual build-up of the flatter, linear portion 
of the Mhh-Msph relation. 

The aim of this work is to use MORGANA to deepen our un¬ 
derstanding of the origin and evolution of the correlation between 
BHs and their hosts, by searching for the key physical processes 
that have an effective impact in shaping the Mhh-Msph relation and 
its proposed break. More specifically, we will explore the impact 
of key physical mechanisms (like disc instabilities, galaxy mergers 
and stellar/AGN feedback) on the observables of interest. MOR¬ 
GANA relies on quite general assumptions on how galaxy forma¬ 
tion takes place inside dark matter haloes, and as such it is rep¬ 
resentative of the whole class of semi- analytic models (SAMs) of 
galaxy formation a nd evolution (see e.g iFontanot et alj2009Ll20l^ : 
iKnebe et al.ll2015| . for comparison with other models in the litera¬ 
ture). 

In the following, we recalibrate the MORGANA model to re¬ 
produce the most recent determinations of the local Mhh-Msph re¬ 
lation and BH mass function, together with the AGN luminosity 
function along the cosmic history. We confirm the predicted break 
in the Mhh-Msph relation, and c heck that it is compatib le with the 
determinations of iGrahanil (1201 2h an d lScottetal.l ( l2013l) . We inves¬ 
tigate the origin of the break, and identify stellar feedback in star¬ 
forming bulges as the mechanism responsible for it. This points to¬ 
wards a co-evolution scenario in which stellar feedback shapes the 
low-mass Mhh-Msph relation in a way that is incompatible with a 
mere co-habitation scenario. 

The structure of this paper is as follows. We first summarise 
the key aspects of the AGN modelling in MORGANA in Sec. 
and detail the different model variants considered in this work in 
Sec. [ 3 ] We then c ompare in Sec |4] model predictions with both the 
IScottet aP (l2013h data on the Mbh-Msph relation, and other phys¬ 
ically linked observational constraints such as the bolometric AGN 
luminosity function and the local BH mass function. We then dis¬ 
cuss the relative contribution of the different physical mechanism 
included in our modeling, and present our conclusions in Sec.|5] 


2 SEMI-ANALYTIC MODEL 

MORGA NA, first presen t ed iniMonaco et al.l ( l2007h and further up¬ 
dated in ILo Faro et al.l (I 2 OO 9 I) . relies on simplified modeling of 
those physical processes that are believed to take place inside dark 
matter haloes and drive the formation of galaxies. These can be 
broadly divided into gravitational processes, such as stellar strip¬ 
ping, secular evolution, galaxy interactions and mergers, and hy- 
drodynamical/thermal processes, such as gas cooling, star forma¬ 
tion, BH accretion, stellar and AGN feedback. The main strength 
of the semi-analytic approach is the extensive use of approximate 
prescriptions, based on theoretical, numerical, or observational re¬ 
sults, to follow complex physical processes in a simplified way. 

We refer t he reader to the original papers and to our compar¬ 
ison work (e.g. Ipontanot et al.ll200^ : IPe Lucia et alj|20l^ l2(}lll : 
iFontanot et al] |2013l) for a complete overview of the model; in 
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the next section we will mainly focus on a single aspect of MOR- 
GA NA, namely the modelling of BH accretion and AGN feedback 
isee lMonaco & FontanoJ200a : lFontanot et al.l2006L for a complete 
overview of this approach). 


2.1 BH accretion 


The m odel for accretion of gas onto the BH follows lOranato et al.l 
( |2004|) . and starts from the assumption that the main bottleneck is 
given by the loss of angular momentum, necessary for the gas to 
flow onto a putative accretion disc. As soon as a galaxy is formed, 
the model assumes that it contains a seed BH of 10® Mq. The first 
step in the loss of angular momentum is connected to the same pro¬ 
cesses (mergers, disc instabilities) that bring gas to the spheroidal 
component; this assumption links BH growth with bulge formation 
events in the history of the model galaxy. Further loss is triggered 
by other physical processes (i.e. turbulence, magnetic fie lds or radi¬ 
ation drag) that are related to star formation in the bulge dUmemural 
l200lh . This results in the building up of a reservoir of low angu¬ 
lar momentum gas with mass Mrs. Its evolution is regulated by 
growth and loss rates (Mj^g and Mj^g respectively). The growth 
rate is taken to be proportional to some power of the bulge star 
formation rate (0b): 

where a = 1 refe rs to the lUmemur i ( I 2001 I) relation used in 
lOranato et al.l ll2004l . while a = 2 corresponds to a model 
where the loss of an gular momentum is due by cloud encounters 
dCattaneo et al.l200^ . /bh is a free parameter, whose value affects 
the normalization of the final Mbh-Msph relation, as it regulates the 
total amount of cold gas flowing into the reservoir during galaxy 
evolution. The gas in the reservoir then accretes on to the BH at a 
rate r egulated by the viscosity of the accretion disc iOranato et al.l 
l2004h : 


where as is the ID velocity dispersion of the bulge and Mbh is 
the BH mass. This term is also capped at the Eddington limit. 

A key role in the feeding of the central BH is played by the 
actual amount of cold gas available in the spheroidal component 
as a consequence of bulge formation events. A comprehensive dis¬ 
cussion of the relative importance of the different physical mecha¬ 
nisms responsible for the buil ding up of spheroidal c omponents in 
galaxies has been presented in lPe Lucia et al.l 1 I 2 OIII) . 

Here we simply recall that in MORGANA galaxy mergers are 
distinguished in major (baryonic merger ratio larger than 0.3) and 
minor mergers. In the former case, the whole stellar and gaseous 
contents of both colliding galaxies are transferred to the spheroidal 
remnant, while in the latter the total baryonic mass of the satellite 
is given to the bulge of the central galaxy. 

Secular evolution is an additional mechanism responsible for 
transporting material to the centre o f galaxies. In MORGANA we use 
the standard fEfstathiou et alj ( Il982h criterion to define the stability 
of disc structures against self gravity. However, whenever a disc is 
found to be unstable, no consensus has been reached yet on how to 
model the resulting instability. Different SAMs have made differ¬ 
ent c hoices, from moving just enough materia l to restore stability 
(e.g.. [Guo et al.ll201ll : lHirschmann et al.l 2012h. to completely de¬ 
stroying the unstab l e disc (e.g.. Bower et'^lTOOd ; lEanidakis et al.l 
l20lA . fMenci et al.l l l20l4 used an updated version of a model by 


iHonkins & Ouata^ ( 1201 ll) . calibrated on results from numerical 
simulations. In all cases, whenever an unstable disc contains cold 
gas, a fraction of it will be made available to accrete onto the BH, 
powering an AGN. In MORGANA, we assume that, when the insta¬ 
bility criterion is met, a fixed fraction (/di = 0.5) of the baryonic 
mass of an unstable disc goes into the bulge component. Cold gas 
flowing into the bulge can then accrete onto the BH as described 
above. In what follows we will discuss the impact of varying /di 
on our final results. 


2.2 AGN-triggered winds 

AGN activity can inject a relevant amount of energy into the 
inter-stellar medium, favouring the triggering of a massive galac¬ 
tic wind that halts the star formation episode ultimately responsible 
for BH accretion. In the AGN feedback model included in MOR¬ 
GANA, the onset of such dramatic events is supposed to be due to 
the interplay between stell ar and AGN feedback, as described in 
[Monaco & Fontaii^ ( l2005h . The triggering of a massive wind ca¬ 
pable of removing all cold gas from the host bulge, is subject to 
three conditions. The first condition is motivated by the theoretical 
expectation that AGN radiation is able to evaporate some 50Mq 
of cold gas for each Mq of accreted mass. Whenever the evapora¬ 
tion rate overtakes the star formation rate, the inter-stellar medium 
is expected to thicken and cause the percolation of all Supemovae 
remnants into a galaxy-wide super-bubble, that can be accelerated 
by radiation pressure from the AGN. The second condition requires 
that energy needed to sweep out all the inter-stellar medium of the 
host bulge is not greater than the available energy from the AGN. 
The third condition is that the BH is accreting in a radiatively effi¬ 
cient w ay, i.e. more tha n one p ercent of its Eddington limit. 

In iFontanot et al.l l l2006h . we showed that AGN-triggered 
winds are a key ingredient to reproduce the evolution of the AGN 
luminosity function and the space density of bright quasars in MOR¬ 
GANA: therefore, in this paper we will consider models including 
this mechanism. In particular, we will refer to the MORGANA ver¬ 
sion labelled “dry winds”, which implies using a high value of /bh 
with winds actually limiting the BH masses. It is worth noting, 
that AGN-triggered winds are not necessary to reproduce the Mbh- 
Msph relation in MORGANA. In a model without winds, however, 
the Mbh-Msph relation is mainly determined by the mechanism 
that regulates the accretion onto the BH, rat her than by feedback- 
based self-regulation. We refer the reader to lFontanot et al.l l l2006h 
for more details, here we just stress that, within MORGANA, the 
AGN-triggered winds have a limited impact on the shape of the 
Mbh-Msph relation (their Fig. 6), although they can affect its red- 
shift evolution (their Fig. 10). 


2.3 Role of stellar feedback 

The stellar feedback prescriptions includ ed in MORGANA follow 
the results from the analytic model of iMonac 3 (l2004l) . which 
postulates different regimes according to the density and vertical 
scalelength of the galactic system and distinguishes “thin” sys- 
tems (i.e discs) from “thick” systems (i.e. bulges and spheroids). In 
IFontanot et al.l ( l2006l) . we showed that stellar feedback has an over¬ 
all modest impact on the evolution of the AGN population, with the 
relevant exception of kinetic feedback in thick systems. 

Kinetic stellar feedback is modelled by estimating the level of 
turbulence in a star-forming bulge, quantified as the gas velocity 
dispersion acg, due to the injection of kinetic energy from Super¬ 
novae and the dissipation of turbulence on a sound crossing time 
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Figure 1. Redshift evolution for the AGN bolometric luminosity function. 
Yellow areas represent an estima te for the bolome tric luminosity function, 
obtained from the Hard X-ray LF iUeda et al.l2014l) . Black solid, dashed red 
and dot-dashed blue lines refer to the predictions of the reference, low-DI 
and high-DI runs respectively. Hatched areas show the l-tr variance associ¬ 
ated with model predictions for the reference model. 


scale. It is easy to show that this quantity scales with the gas con¬ 
sumption time-scale of gas f* to the power —1/3, so the velocity 
dispersion of gas in star-forming bulges can be modeled as follows: 


(Teg 


Co 



- 1/3 

km/s 


(3) 


where the paramete^^ Co regulates the efficiency at which Super¬ 
novae drive turbulence. Assuming that cold gas clouds have a 
Maxwellian distribution of velocities, cold gas is then ejected into 
the halo at a rate proportional to the probability that the velocity of 
a cloud overtakes the escape velocity from the bulge. So, the com¬ 
bination of a significant and a lower f* leads to high values of 
(Teg in bulges, driving massive outflows in small bulges. 

In thin systems, where Supernova energy can be easily blown 
out of the sys t em, th e value of (o is expected to be smaller. In 
iFontanot et al.l ( l2006h we showed that kinetic feedback mechanism 
in star-forming bulges is a viable mechanism to explain the down¬ 
sizing behaviour of the AGN population, as high velocity disper¬ 
sions of cold gas leads to massive removal of cold gas and to the 
suppression of low-luminosity AGN at moderate-to-high redshifts. 


3 RUNS 

The reference AGN feedback model by iFontanot'et ^ 1 I 2 OO 6 I 1 
adopted here has two main parameters that regulate the accretion 
rate of the cold gas reservoir around the BH, namely /o, which reg¬ 
ulates the amount of stellar kinetic feedback, and /bh which tunes 


2 


This parameter was named (tq in 


IFontanot et alj lioo^i . 


Table 1. AGN feedback parameters for the main runs discussed in this work. 


Model 

/di 

Co [Km/s] 

/bh 

Reference 

0.5 

30 

0.05 

Low-DI 

0.0 

20 

0.02 

High-DI 

1.0 

20 

0.003 


the normalization of the Mbh-M^ph relation. In addition, the over¬ 
all radiative efficiency rj, which is the emitted luminosity in units 
of the rest-mass energy, must be recalibrated when using a higher 
normalization for the Mbh-AFsph relation. This quantity is in fact 
broadly proportional to the ratio between the integrated emissiv- 
ity of A GN over tim e and luminosity and the local mass density 
of BHs (ISoltanlll982h . We thus decrease it from the original value 
of Tj = 0.1 adopte d by IFontanot et al.l ( l2006h to 7 ; = 0.06 (see 
IShankar et al.ll2013h . and keep it always fixed to this value in what 
follows. 

All the runs presented in this work have been performed 
on the same cubic box (100 Mpe side) obtained using the code 
PINOCCHIO jMonaco et alj|2003.1201 3h with N = 1000® parti¬ 
cles and assuming a 2006 concordance cosmology with parameters 
flo = 0.24, Qa = 0.76, h = 0.72, as = 0.8, Usp = 0.96. 


As the original MORGANA runs were calibrated agains t 
the observed Mhh-Msph relation from iMarconi & Hunt! (l2003b . 
we have recalib r ated t he model to fit the relation proposed by 
iGraham & Scotj l l2015h {reference run), and to reproduce the evo¬ 
lution of the AGN luminosity function. Regarding the Mhh-Msph 
relation, we decided to calibrate our models to reproduce the most 
reliable BH mass estimates, in the bulge mass range of M* > 
10®® M 0 . The best-fit parameter values are reported in Table [T] 
Models tend to have marginally steeper slopes in the high mass 
end of the Mhh-Msph relation, so in general the relation is under¬ 
estimated at the bending point. It would be easy to absorb this dis¬ 
crepancy using a different calibration choice; however, we consider 
this conservative choice sufficient for the scope of this paper. 

We then define some model variants, to explore the sensitiv¬ 
ity of our results to model assumptions. We first explore the im¬ 
pact of our choice for the modelling of bar instabilities, by defining 
two additional runs. In the first one, we switch off disc instabilities 
(/di = 0, no-DI run), while in the other, we assume that the un¬ 
stable disc is completely destroyed and the disc instability results 
into a purely elliptical galaxy (/di = 1 high-DI run ) . The se runs 
have also been recalibrated to fit the lGraham & Scol^ ( 1201 Sl) Mbh- 
Msph relation. Best-fit values for the relevant parameters have been 
collected in Table[T] 

The AGN bolometric luminosity function used for calibration 
is reported in Fig. [T] Here the reference run is shown as a black 
solid line, while dashed red and dot-dashed blue lines refer to the 
no-DI and high-DI runs respectively. The observational estimate of 
the bolometric luminosit y function (yellow area) is obtained from 
the Hard X-ra y LF from Ueda et all (I20l4 yellow shaded area), 
using both the lMarconi et al. U2004t) bolometric correc tion and an 
estima te for the fraction of Compton-thick objects by lUeda et al.l 
( l2014t) . All three runs broadly reproduce the evolution of the bolo¬ 
metric AGN-LF over wide ranges of redshifts and luminosities, al¬ 
though the high-DI model tends to overproduce AGNs, especially 
faint ones, at high redshift. It is worth noting that the new cali¬ 
brations assume relative l y low values for /o with respect to those 
quoted in iFontanot et al.ld2006h . This implies that the impact of this 
parameter to regulate the downsizing behaviour of the AGN popu- 
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Figur e 2. Mth-A^sph relation at z = 0, for differe nt models considered in t his work. Red lines correspond to the best fits to the observa tions as inlScotLsLal 


20U ), while the green line refers to the best fit as in iKormendv & h 3 j2013h : open circles and squares show data from IScott et alj i2013h and iGraham & Sc^ 


201511 . respectively. Solid black lines show the median relation for model galaxies, while the grey area represents the 16% and 84% percentiles. 


lation is reduced and also runs with C^o = 0 (see below) provide a 
fair agreement with the data. 

For each of the three disc instability options we perform addi¬ 
tional MORGANA runs by switching off the stellar kinetic feedback 
((^0 = 0). These model variants have not been recalibrated with 
respect to the parameters used in the > 0 runs. A further test 
has been performed by running the model without using the cold, 
low angular momentum gas reservoir, i.e. assuming that a (Edding¬ 
ton limited) fraction Mj[g of the cold gas in the bulge is directly 
accreted by the BIT. Finally, we have studied the role played by 
galaxy mergers on the Mbh-ATsph relation. 

In MORGANA, the merger times for satellite galaxies^ (fmrg) is 
co mputed using an upd ated version of the fitting formulae pro vided 
by iTaffoni et al J ( l2003ll . As shown in lPe Lucia et al.l l l2oi(]|) . over 
the relevant range of mass ratios, these fmrg are typically shorter 
than those estimated from other authors, using different and more 
recent numerical simulations. We then present two model variants 
obtained using different fmrg definitions, namely the fitting formu- 

® Galaxy merger times are computed whenever a galaxy becomes a satel¬ 
lite, following the merger of its parent DM halo with a larger one. 


lae proposed by iBovlan-Kolchin et al.l 1 I 2 OO 8 L “longer fmrg”) and 
the more extreme case of instantaneous mergers (fmrg = 0). Also 
in these three additional versions we kept the other parameters as 
in the ^0 > 0 runs. 

It is also worth stressing, that in all the models discussed in 
this work we do change only the subset of parameters directly con¬ 
nected to AGN mod elling, while k e eping all other parameters to 
the values defined in ILo Faro et alj ( 1200^ . with the exception of 
the recomputation of the size of starburstO 

This implies that, in general, the agreement of model predic¬ 
tions with the calibration set (i.e. local stellar mass function, cosmic 

^ The computa tion of gas consumption time-scale t*, based on the 
lKennicut3 il998ll relation and thus on gas surface density, requires knowl¬ 
edge of the size of the star-forming region. Thi s quantity was set equ al to the 
bulge size in the original MORGANA, while in lLo Faro et in ( 2 ^ a more 
refined modeling was adopted, based on the idea that the level of turbulence 
should equate the velocity given by the bulge rotation curve computed at 
the starburst size. Nonetheless, this change caused a drastic steepening of 
the Mbh-^sph relation, so we switched it off to achieve a good modelling 
of AGNs. The impact on other observables is modest, and can be absorbed 
by suitable retuning of other parameters. 
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Figure 3. Mbh-Afsph relation at z = 0, for different models considered in this work. Lines, colours and shading as in Fig.|2] 


star formation rates and morphological mix) may not necessarily be 
optimal. We, however, explicitly tested that the deviations from the 
reference models for non-AGN galaxy properties are small. For in¬ 
stance, stellar mass functions are relatively weakly affected by the 
exact AGN modelling adopted, once the model is calibrated to re¬ 
produce the Mhh-Maph relation and the AGN luminosity function, 
while the morphological mix is more sensitive to the (^o parameter 
and to disc instabilities. We checked that our reference run produces 
a plausible morphological mix. 


4 RESULTS 


The z — 0 Mhh-Maph relation for the runs defined in the previ¬ 
ous section are shown in Fig.[2 In all panels, the black solid line 
is the median relation for the model galaxy sample, while the grey 
area represents the 16th percent and 84th percent percentiles of the 
distribution. The central panel in the upper row displays the Mbh- 
Msph relation in our reference run, compared to the observational 
relations proposed by Kormendv & Hoi ( 1201 3l . green solid line) and 
iGraham & Scotj ( 120151 red solid line). A steepening of the Mbh- 
Msph relation is evident for M* 10 ^° Mq gala xies, resulting in 
a brea k of the relation similar to the results of iGraham & Scota 
( l2015h (red lines represent their best fit to power-law-Sersic and 
core-Sersic samples). Moreover, the scatter in the mode l is broadly 
compatible with that in the data. The lScott et alj ( 1201 3h datapoints 
suggest a larger scatter than the predicted one. Note, however, that 
our model does not include observational errors on BH and stel¬ 
lar ma sses. Also, the intrinsic scatter reported bv iKormendv & ll3 
( l2013h for the subsample inclusive of only the most accurate BH 
mass estimates, amounts to ~ 0.3 dex, even smaller than what pre¬ 
dicted by our models for the total sample. 

We now discuss the implications for the Mhh-Maph relation 
when varying one or more of the input assumptions in our refer¬ 
ence model. We first consider the effect of changing our modelling 



Log(MB„ [MJ) 

Figure 4. BH m ass function at z = 0. Data from lShankail l2013i grey area) 
IShankail i2009l green hatched area). Linetypes and colours as in Fig.[T] 


for disc instabilities (Fig. [2] upper row). We remind that the no- 
DI and high-DI runs require recalibration of the i^o and /bh pa¬ 
rameters to reproduce the observed Mhh-Maph relation, due to the 
different predicted amounts of cold gas available for BH accretion 
and AGN feedback in the models. For all the runs we find a clear 
break in the Mhh-Maph relation, thus concluding that its presence 
does not depend on the exact modelling of disk instabilities. The 
scatter, compared to the reference run, decreases in the no-DI case. 
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Figure 5. Redshift evolution of the overall Afbh-A^sph relation. Line types 
and colours as in Fig.[T] See text for more details on the definition of the F 
parameter. 


but remains constant in t he other case . This is at variance with the 
results presented, i.e. bv iMenci et al.l (l2014h . We will deepen into 
this point in Section|5] 

We then turn to the predicted BH mas s function (Fig. |4), and 
compare it with the empirical estimates by IShank^ ( l2013h . based 
on th e convolution between the galaxy velocity dispersion function 
from lBemardi^t_j! l201(]|), a nd the BH mass-velocity dispersion 
relation by iMcConnell & m 3 ( l2013h . which is the most appropri¬ 
ate for our renormalized-high BH-stellar mass rela t ion. F or com¬ 
pleteness, we also report the estimates by IShan ( l2009ll . which 
were anyway based on previous (lower) normalization of the BH 
scaling relations. All our runs show a reasonable agreement with 
the most recent estimate at the high-mass end. At low BH masses 
the models with no or moderate DI tend to underpredict the mass 
function, which was anyway estimated assuming a strictly linear 
power-law scaling relations. The model with strong DI provides 
a better fit, at the cost of an overprediction at large massefl We 
caution that the comparison with the local BHMFs is simply pro¬ 
vided as a broad consistency check. As discusse d by several authors 
(see e.g. lTundo et al.l2007l : fshankar et alj2012h . these estimates are 
still affected by systematics on the exact normalization of the input 
scaling relations and on the actual contributions of low-mass black 
holes, which can be connected either to possible breaks or other 
deviations (like, e.g., those proposed for the “pseudo”-bulge pop¬ 
ulation) from the scaling relations defined by more massive black 
holes. 

The results shown in Fig.[T]and|4]confirm our previous conclu- 


® This trend cannot be predicted by the results of Fig.|^ because it is driven 
by the larger mass in bulges that is obtained in this model. A full recalibra¬ 
tion of the model, including also galaxy observables, may absorb this dif¬ 
ference, although such a strong role of disc instabilities may easily result in 
an overproduction of bulges. 


sion that models with different strength of disk instability provide 
very similar descriptions of the evolution of the AGN population. 

We then use our additional model variants to understand which 
physical mechanism implemented in MORGANA is the main driver 
for the break in the theoretical Mhh-Msph relation. In the lower 
row of Fig. [2 we show MORGANA runs where we switched off 
kinetic stellar feedback in bulges (i.e. i^q = 0). All the resulting 
relations are power laws over the whole Msph range, and follow 
a relation with a slope c ompatible to, or slightly steeper than, the 
iKormendv & Hoi ( 1201 3h best-fit relation (green line), with no ap¬ 
parent bend at the low-mass end. We thus conclude that the treat¬ 
ment of stellar feedback is the main responsible for the different 
behaviours at the high- and low-mass end of the Mhh-Msph rela¬ 
tion. 

This is confirmed in Fig. [3 where we report the Mhh-Msph re¬ 
lation predicted by the variants of the reference model without the 
reservoir of low angular momentum gas, and with varied merging 
times. We recall that these additional variants do use the same cali¬ 
bration as the reference run, so the normalization of the Mhh-Msph 
relation is not guaranteed to be reproduced. All relations show the 
same steepening at low masses. Moreover, the typical Msph mass 
scale corresponding to the resulting break does not change among 
the various runs (with a possible exception of the zero merger time 
run), confirming fhat sfellar feedback is sefting fhe mass scale of fhe 
change of slope. Changing fhe merger timescale in MORGANA has 
only a limited effect on the shape and normalization of the Mbh- 
Msph relation, and we checked that this is independent of the as¬ 
sumed disc instability model. 

We finally address fhe redshifl evolution of the Mhh-Msph 
relation. We checked that all models predict the Mhh-Msph relation 
to be already in place at high redshift, with almost no evolution in 
the shape with respect to the local relation. It is worth stressing that 
the actual detectability of a “bended” relation at higher redshifts is 
related to the abundance of the massive galaxy population. 

We find some evolution in the normalization of the Mhh-Msph 
relation. We quantify this evolutionary trend by means of the r(z) 
parameter, defined as: 


r(.) = 

■^sph 


Mbh(z) 
Mbh(2 = 0) 


(Msph) 


(4) 


In Fig. |5] we show the average evolution of F, computed with 
Msph > 10® Mq model galaxies, using the corresponding median 
relations shown in Fig. [2] as 2 = 0 reference. The evolution up 
to 2 = 2 is similar for all models, but disc instabilities lead to a 
flattening of F at high redshift. 


5 DISCUSSION & CONCLUSIONS 


Recent re-calibrations of the locally measured Mhh-Msph relation 
suggest its normali zation to be a factor of ~ 2 higher than previ - 
ously inferred (e.g., Gr^imi_^Scotj20]A ; Kormendv & Hol2013l) . 
Some authors dScott et al.l2013l : lGraham & Scottl2015h have found 
a strong, quadratic steepening of the Mhh-Msph relation below 
Mbh 2 X 10® Mq. The MORGANA mo del made a predictio n 
of such a steepening at low masses in 2006 dFontanot et al.ll200^ . 
While consensus on this break in the relation is still to be achieved, 
in this paper we have further investigated the possible physical 
causes of such a steepening. To this aim, after recalibrating the pre¬ 
dicted Mhh-Msph relation to take account of the higher normal¬ 
ization, we explored the impact of varying the strength of the disc 
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Figure 6. Additional black hole - spheroid scaling relations for at 2 = 0, for the different models considered in this work. Left panel: M|-,h-A^sph relation 
relation (data as in Fig.|2j. Middle panel: Mbh-^dm relation. Right panel: Mgph-MDM relation. In all panels, solid lines refer to the median relations in 
reference model (shaded areas representing the the 5% and 95% percentiles), while the dashed line to the corresponding Co = 0 runs. 


instabilities, galaxy mergers and stellar feedback in star-forming 
spheroids. 

We explicitly tested that the strength of the disc instability, i.e., 
the amount of baryonic material transferred from the unstable disc 
to the central spheroidal component, does not affect the shape of the 
Mhh-Maph relation, although it may change its normalization. We 
also showed that neither the assumed modelling of a gas reservoir 
around the BH, nor the timing of galaxy mergers have a relevant 
effect on the shape of the Mbh-A/gph relation in MORGANA, while 
they have an influence on its normalization and scatter. 

In more detail, we find that the scatter in the Mhh-Msph rela¬ 
tion, as predicted by MORGANA does not depend on the modeling 
of disc instabilities, i.e., on the relative importance of this physical 
mechanism in feeding the central BH. This re sult may be recon - 
ciled with the results from other SAMs (see e.g lMenci et alj2oT^ . 
by considering the role of the gas reservoir in our model. As a mat¬ 
ter of fact, comparing the runs assuming direct (Eddington limited) 
accretion onto the central BH, i.e., switching off the modeling of 
the gas reservoir, we confirm a decrease (of a factor about 1.5) in 
the scatter at the high-mass end of the relation moving from the 
reference to the high-DI runs. This is compatible with the idea that 
the delayed accretion from the reservoir (with respect to the time 
when cold gas looses most of its angular momentum) adds another 
time scale to the mechanism of BH accretion that couples to disc 
instabilities. This additional time scale dominates the scatter in the 
Mhh-Maph relation. 

The physical process responsible for the break in the Mbh- 
Maph relation in MORGANA is stellar feedback in star-forming 
bulges. Its role can be illustrated as follows. Both bulges and 
BHs acquire stars from merging and local star formation/a ccretion. 
While most stars are found in bulges at 2 ; = 0 (see e.g. iGadottil 
l2009l) . star formation mostly takes place in the disc components of 
galaxies, with a minor contrihution from starhursts. This is true in 


MORGANA jMonaco et al.ll200% . but observations suggest a sim¬ 
ilar trend at 2 ~ 2, with normal star-forming galaxies being pre¬ 
dominantly rotating discs jporster Schreiber et alj|20()^ and star- 
bursts lying above the main sequence of starforming g alaxies, con¬ 
tribut i ng only ^10 per cent t o the total SFR densi ty l lDaddi et al.l 
I 2 OIOI : iRodighiero et al.l 1201 ll : ISargent et al.l l2012h . As a conse¬ 
quence, most stars in bulges were formed in discs and were car¬ 
ried into the spheroidal component by mergers and disc instabili¬ 
ties. Conversely, m ost mass in the BHs is brought by accretion and 
not by mergers (e.g. ISalucci et al.ll999l:lMarconi & IIunll20()l . As 
long as BH accretion is related to star formation in bulges, a selec¬ 
tive suppression of this star formation in small bulges will decrease 
small BH masses but will not influence much the stellar mass of 
bulges. 

This interpretation of the b reak Mhh-Maph relation differs 
from the ICraham & Scotll (l2015h . and we propose it as an alter¬ 
native explanation of the steepening of the Mhh-Maph relation 
at small masses. In particular, our model does not produce two 
clearly different regimes as a function of Maph- Star formatioifl 
BH accretion and mergers (both dry and wet) happen at all mass 
scales, although the relative importance of each physical process 
depends on the tota l mass of the host galaxy dPe Lucia et alj20TII : 
IShankar et al.l201^ . We also showed that the the shape of the Mbh- 
Msph relation depends weakly on the choice of merger timescales. 
In particular, a significant shortening of the merger timescales (i.e., 
assuming instantaneous galaxy mergers) increases the mass of the 
most massive BHs by a factor of up to ~ 2 and changes the position 


® In detail, star formation is still appreciable for massive galaxies in 
MORGANA, due to the inefficient quenching of cooling flow in massive 
haloes for an AGN m odelling, which does not include hot gas accretion 
dPontanoT et alj200% . 
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of the break, but the overall Mbh-A/gph relation still shows a bend 
at the low-mass end. 

We also stress that the stellar feedback mechanism we pro¬ 
pose is an example of a possible self-regulation of BHs and bulges, 
where effective stellar feedback leads to starvation of small BHs. 
In this sense, a steepening of the Mhh-Maph relation can be seen as 
an argument in favour of the “co-evolution” scenario as opposed to 
the “co-habitation” one, where a linear relation driven by mergers 
is the most natural expectation. This self-regulation is not due to 
AGN feedback, although, as commented in Section[2 the final BH 
mass is self-regulated by A GN feedback. Indeed, of the three mod¬ 
els presented in the original lpontanot et alj j2006h paper, two were 
based on a self regulation of BH masses due to AGN feedback, but 
all three showed a break in the Mhh-Msph relation. 

Previous work already showed that non-trivial AGN feedback 
schemes, whose efficiency varies as a function of halo mass, can 
creat e a break in scaling relations, including the Mbh-ATgph rela¬ 
tion Jcirasuolo et "^bOOSl : IShank^ et al.ll20 06h. though p ossibly 
milder than the one emphasised bv iGraham & Scod l l2015l) . 

It is at this point natural to ask whether a break in the BH- 
stellar mass relation corresponds to a break in other scaling rela¬ 
tions, especially between M bh and velocity d i spersi on a, usually 
measured to be even tighter. iGraham & Scotj ( 1201 Sll have shown 
that no apparent break is evident in their data when plotting BH 
mass versus velocity di spersion, while the Su pernova feedback 
model implemented by ICirasuolo et al.l ( l2005h predicted a clear 
break also in this relation. An estimate of the central velocity dis¬ 
persions of bulges from MORGANA would require to make several 
additional assumptions on the exact light and mass profiles of the 
stellar, gas, and dark matter components. Instead, we choose to con¬ 
sider the predicted relation between Mbh and halo mass Mum- 
this is shown in Fig. [h] where the left panel is a replication of the 
reference model relation in Fig. [^(continuous line) and, overplot¬ 
ted, the median relation for the model with (^o = 0 (dashed line). 
The other two panels of this figure show the relation between BH 
mass and halo mass (middle), and between halo mass and stellar 
mass of the spheroidal component (right). A break in the Mbh- 
Mdm relation is evident in our reference model, while it is nearly 
absent in the model without kinetic feedback (dashed lines). It is 
interesting to notice that the shape of the Msph-MoM relation is 
almost unaffected by this mechanisms, confirming that this mod¬ 
elling of feedback in star-forming bulges affects BH masses much 
more than the spheroidal component of the galaxy. 

It is worth stressing that it is not straightforward to infer con¬ 
clusions on the shape of the Mbh- o' relation from the results shown 
in Fig. While we expect Mdm to correlate with virial velocity, 
the relation between the latter quantity and velocity dispersion is 
still uncertain for large galaxy samples. S everal groups (see e.g. 
iFerrarese & Merritll20od ' lBaes et alj2003l) suggested that a corre¬ 
lates with the circul ar velocity as measured at the outermost radii, 
while other authors i IHoI| 2007|) showed that such a correlation has a 
significant scatter and varies systematically as a function of galaxy 
properties. Therefore, any further conclusion on the Mbh-ct rela¬ 
tion critically depends on the assumed relation between velocity 
dispersion and circular velocity, e.g. if we assume a linear correla¬ 
tion, we then expect a break also in the this relation. 

Of course, our modelling of stellar and AGN feedback still 
represents an idealized approach and we argue tha t observations 
of spatially resolved galaxies (like e.g. CALIFjAQ, I Sanchez et al.l 


^ Calar Alto Legacy Integral Field Area Survey 


l2012h will provide a better comprehension of the interplay between 
these two fundamental physical mechanisms. 
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